Photon polarization and geometric phase in general relativity 
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Rotation of polarization in an external gravitational field is one of the effects of general relativity 
that can serve as a basis for its precision tests. A careful analysis of reference frames is crucial for 
a proper evaluation of this effect. We introduce an operationally-motivated local reference frame 
that allows for a particularly simple description. We present a solution of null geodesies in Kerr 
space-time that is organized around a new expansion parameter, allowing a better control of the 
series, and use it to calculate the resulting polarization rotation. While this rotation depends on the 
reference-frame convention, we demonstrate a gauge-independent geometric phase for closed paths 
in general space-times. 

I. INTRODUCTION 

Electromagnetic waves — visible light and other bands of the spectrum — are our prime source of information about 



QfMcctromagncuc waves — visiDie ngnt ana otner Danas or tne spectrum — are our prime source oi nuormaxion apout 
the Universe Since measurements of the light deflection near the Sun [H were made in 1919 wave propagation 

OO ' is used to test general relativity (GR). Two of the "classical tests" of GR, light deflection and time delay, can be 

| understood in terms of geometric optics ■ The first post-eikonal approximation [3| allows to track the evolution of 

^ ■ electric and magnetic fields along the light rays and thus discuss polarization. 

' In this approximation we speak about photons with a null four-momentum k and a transversal four-vector polar- 

O" 1 ization f. Both vectors are parallel-transported along the trajectory, which is a null geodesic 0, @: 

feb: k-k=o, v k k = o, (i) 

k-f = 0, V k f = 0, (2) 

<N ■ 

■ where Vk is a covariant derivative along k. 

| In the last decades polarization has been yielding important astrophysical and cosmological data. Cosmic microwave 
t"**- ■ background Q, blazar flares @,[l^, astrophysical jets [tlGil, and searches for dark matter [ll[ are just several examples 
ON ■ where polarization conveys crucial information. 

, Photons are commonly used as physical carriers of quantum information. The abstract unit of quantum information 
is a qubit (a quantum bit) [lH, and two linearly independent polarizations encode the two basis states of a qubit. One 
| of the branches of quantum information is quantum metrology, which aims to improve precision measurements by 
using explicitly quantum effects, such as entanglement (l3| . Relativistic properties of the information carriers [Til [l5| 
become important whe n q uantum technology is used in precision tests of relativity (l6j or quantum information 
processing on the orbit p7| . 

Gravity causes polarization to rotate. This effect is known as a gravimagnctic/Faraday/Rytov-Skrotskii rotation 
(l8l - l2i| . Helicity is invariant under rotations, but states of a definite helicity acquire phases e ±lAx . Depending on the 
context we refer to Ay either as a polarization rotation or as a phase. Once evaluated this phase can be encapsulated 
as a quantum gate [12j and incorporated into (quantum) communication protocols or metrology tasks. 

In the Schwarzschild space-time, as well as at the leading order of the post-Newtonian approximation, this phase is 
known to be zero [HI, [IH, [H, [23| . It is higher-order gravitational moments that are held responsible for the rotation 
of the polarization plane. In particular, the GR effects were shown to dramatically alter polarization of the X-ray 
radiation that is coming from the accretion disc of the (then presumed) black hole in Cyg X-l [24| . Numerous 
analytical and numerical studies of trajectories and polarizations in different models and astrophysical regimes were 
performed ( p~9j - [2ll [25M29| . and references therein). The scenarios included fast- moving gravitating bodies, influence 
of gravitational lenses and propagation through gravitational waves. 

The results are often contradict ory . Some of the contradictions result from genuine differences in superfluously 
similar physical situations [20l l28l . [29{ ■ On the other hand, polarization rotation is operationally meaningful only if 
the evolving polarization vector is compared with some standard polarization basis (two linear polarizations, right- 
and left-circular polarizations, etc.) at each point along the ray (30. |32|. 

Setting up and aligning detectors requires alignment of local reference frames. Much work has been done recently 
on the role of reference frames in communications, especially in the context of quantum information |3lj . Partial 
knowledge of reference frames can lead to loss of communication capacity, and mistakes in identifying the information 
content of a physical system. The lack of definition for polarization standards and an ad hoc introduction of the 
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angle adjustments is one reason for the variety of quoted values for the phase A\- Even the consensual result of a 
zero phase in the Schwarzschild space-time should be qualified. Without specifying the appropriate reference frame 
it is either meaningless or wrong. 

In this article we investigate the role of local reference frames in defining A\- The result is obviously gauge- 
dependent, and additional considerations should be used to fix the gauge. Communicating reference frames is some- 
times a difficult procedure which may involve a high communication cost. On a curved background of GR it may 
require some knowledge of the metric at each point along the trajectory. Using the method presented by us in [32J we 
build a local reference frame in stationary space-times, and then fix the standard polarizations in a single construction 
using what we call the Newton gauge. Our construction does not require communication between the parties, gives a 
precise meaning to the idea that there is no polarization rotation in the Schwarzschild space-time and reproduces the 
absence of phase in the Minkowski space-time. From an operational point of view it allows us to set up detectors at 
any point in space based only on the local properties at that point. We illustrate the use of this gauge by studying 
polarization rotation in the Kerr space-time, and derive an explicit expression for Ax in the scattering scenario that 
is described below. The calculations are based on the results of @, [U [2f| with the added value of the Newton 
gauge. Careful bookkeeping is required with respect to the various coordinate systems used in different parts of the 
calculation, as well as different orders of the series expansion. To simplify the latter, we introduce a new expansion 
parameter. 

Closed paths result in a gauge-invariant gravity-induced phase. We discuss them in Section IV. The rest of the 
paper is organized as follows: First we discuss null trajectories in the Kerr space-time. In Section III we review 
the Wigner's construction of setting standard polarizations in a local frame, set-up the Newton gauge and explicitly 
calculate the rotation in the scattering scenario, where the light is emitted and observed far from the gravitating 
body (ri,f2 S> M, a), but can pass close to it. Finally, we discuss gauge-dependent and gauge- invariant aspects of 
our results from mathematical and operational points of view (Section IV). Summary of the important facts about 
the Kerr space-time, as well as detailed calculations and special cases are presented in the appendices. 

We use h H — h signature, set G = c = 1 and use Einstein's summation convention in all dimensions. Three- 
dimensional vectors are written in boldface and the unit vectors are distinguished by carets, such as b. Local tetrad 
components are written with carets on indices, such as and the four- vector itself as k = k^e^, where are 
vectors of a local orthonormal tetrad. 



II. NULL GEODESICS 
A. Null geodesies in the Kerr space-time 



The Kerr metric in the Boyer-Lyndquist coordinates is given by 



2Mr\ , 2 p 2 , , , „ /, , 2Mra 2 2 \ ,„ , , 4Mar 
ds z = - ( 1 — ) dt 2 + t^dr 2 + p 2 d6 2 + ir 2 + a 2 + ^ sin 2 6» ) sin 2 (9d0 2 — shSOdtify, 



where M is the mass and a = J/M the angular momentum per unit mass of the gravitating body. We use the 
standard notations 

p 2 = r 2 + a 2 cos 2 6, A = r 2 - 2Mr + a 2 . (4) 

Thanks to the three conserved quantities — the energy E = — fc , the z-component of the angular momentum L 
and the constant r\ 0, |35[ — the geodesic equations in Kerr space-time are integrable in quadratures. When dealing 
with photons it is convenient to set the energy to unity and to re-scale other quantities [6[. Explicit form and the 
asymptotic expansions of D = L/E and r\ are given in Appendix [XJ Fixing the energy and using the null vector 
condition leaves us with only two components of the four-momentum k. In specifying the initial data we usually take 
k\ and kf as independent. 

We discuss polarization from the point of view of static observers that are at rest in the "absolute" space t = const. 
At every point (outside the ergosphere, if the model represents a black hole) we introduce a chronometric orthonormal 
frame [5|, |3J] (Appendix |A"| . which will be used to express the initial conditions and observed quantities. 

The Hamilton- Jacobi equation for a null geodesic separates and the trajectories can be deduced from it Q. We 
label their initial and final points as (ri,#i,0i) and (r 2 , #2 1 <f>2 ) 1 respectively. The gauge convention that we adapt 
(Sec. MI B[) makes the knowledge of A(j> redundant for discussions of polarization rotation in the scattering scenario 
(For the calculation of <fi see Appendix [X| . In particular we are interested in the integrals 

f r dr f e dd 
31 : = / 7== / 7=' 5 

J ±Vr J ±Ve 
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where 

R = r 4 + (a 2 -D 2 - r,)r 2 + 2M(r) + (D - a) 2 )r - a 2 ^, (6) 
9 = r/ + a 2 cos 2 6 - D 2 cot 2 9. (7) 

Null geodesies in Kerr space-time are classified according to the sign of the constant r\. The scattering scenario 
corresponds to 77 > 0, and typically 77 ~ r\ . Following [25| we also assume that the angle 9 reaches its extremal value 
(either maximum or minimum) only once, 9\ — > # max / m i n — > 9 2 . Series expansions of the above integrals are most 
conveniently written with the help of a new constant 

A 2 :=D 2 + ri. (8) 

We use the parametrization D = A cos a, and 77 = A 2 sin 2 a, with < a < it. In this notation the minimal coordinate 
distance from the center is 

3M 2 2aMcosa a 2 cos 2 a „ 
r min = A-M- — + __ + o(A- 2 ). (9) 

We perform the integration by the methods of 0, HB| (Appendix B). In the scattering scenario 7 , 1 ,r 2 — > 00, while 
the constants of motion are kept finite. Expansion in terms of A gives 

51 = 1 q2 (ipi +1P2 - ° 4 ^ Q (3(^i + -02) + sinipt cos^i + sin^ 2 cos ^2)) + C(A~ 4 )^ . (10) 

Where cost/j = cos0//_t + and /x + is defined in the appendix. 

To separate conceptual issues from the computational details we make two simplifying assumptions. First, the 
initial and final points are taken to lie in the asymptotically flat regions, 7^ — > 00, while 77 and D are kept finite. 
Second, in these regions the polar angle 9 (nearly) reaches its asymptotic values 9 lri and 9 ou t, respectively. 

The asymptotic form of momentum can be deduced from the equations of motion Q. In the chronometric tetrad 
basis the outgoing momentum tends to 

k^(l,l,-^-^-), (11) 
V r r sm 9 out J 



where /3 ou t := +-\/O(0out) and s = ±1. This expression is true in any scenario where an outgoing photon reaches the 
asymptotically flat region. 

The asymptotic expression for the incoming momentum is similar, 

(12) 



where /3j n := Iy9(^), and s = ±1 correspond to 9 ln — > # m ax/#min — > 9 ou t, respectively. 
Equating the radial and angular expressions for 51 we obtain 

<U = - - ia ± jM ^a- cos^ + G(A _ 2 
A sm 9i n 

with the higher order terms and special cases described in Appendix [B] The plus sign corresponds to the trajectory 
in which 9 first decreases with r and reaches # m i n before increasing to 9 ou t- 



B. Null geodesies in 1+3 formalism 

In stationary space-times a tetrad of a static observer is naturally related to Landau-Lifshitz 1+3 formalism (33| . 
Static observers follow the congruence of time-like Killing vectors that defines a projection from the space-time 
manifold M. onto a three-dimensional space S3, 7r : Ai — > £3. 

In practice it is performed by dropping the time-like coordinate of an event, and vectors are projected by a push- 
forward map 7r*k = k in the same way. Contravariant vector components satisfy 



(k) m = (7r*k) m = (k) m = dx m /d\ m = 1,2,3. 



(14) 
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where A is the affine parameter. 

The metric g on M. can be written in terms of a three-dimensional scalar h, a vector g, and a metric 7 on S 3 as 

ds 2 = -h(dx° - g m dx m f + dl 2 , (15) 

where the the three-dimensional distance is given by dl 2 = r y mn dx m dx n . The metric components are 

lrnn=(g mn - § ^A, (16) 



and 

h = -g(X)> Sm = -Eom/Soo- ( 17 ) 

The inner product of three- vectors will always refer to this metric, k f = J mn k m f n - Vector products and differential 
operators are defined as appropriate dual vectors [33$. Finally, the spatial projection of a null geodesic has a length I 
that is related to the affine parameter A as 



dl\ 2 k 2 
dX 



(18) 



For a static observer the three spatial basis vectors of the local orthonormal tetrad are projected into an orthonormal 
triad, 7r*e( m ) = e( m ), e( TO )-e(„) = S mn . We adapt a gauge in which polarization is orthogonal to the observer's 
four- velocity, u-f = 0. In a three-dimensional form this condition reads as kf = 0. 

In stationary space-times the evolution equations Eq. (TTJ), ^ can be reduced into a convenient three-dimensional 
form (20j. Using the relations between four- and three-dimensional covariant derivatives, V M and D m , respectively, 
the propagation equations (TT]) and ^ are brought to a three-dimensional form [ljj [2(| , 

Dk _ , E k, Df _ ~ 

-^nxk + ^k, I = nxf. (19) 

The angular velocity of rotation f2 is 

ft = 2uj - (wk)k- E s x k, (20) 
where g and h play the roles of a vector potential of a gravimagnetic field B 9 , 

u} = -ifc curlg = -ifc B 9 , (21) 



and a scalar potential of a gravielectric field, 



respectively [3£ 



III. LOCAL REFERENCE FRAMES AND POLARIZATION ROTATION 

To discuss how photon's polarization rotates it is necessary to define the standard polarization directions along 
its trajectory. Only after physically defining two standard polarizations it is meaningful to talk about polarization 
rotation, and only in a particular gauge the Schwarzschild space-time induces a zero phase on an open trajectory. 

A. Wigner phase — polarization convention 

Transversality of electromagnetic waves makes the choice of two standard polarization directions momentum- 
dependent. On a curved background it also depends on the location. Wigner's construction of the massless rep- 
resentation of the Poincare group [HI, H3] is the basis for classification of states in quantum field theory. We use it at 
every space-time point to produce standard polarization vectors. 
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The construction — part of the induced representation of the Poincare group |37j — consists of a choice of a 
standard reference momentum kg and two polarizations, a standard Lorentz transformation L(k) that takes kg to 
an arbitrary momentum k, and a decomposition of an arbitrary Lorentz transformation A in terms of the standard 
transformations and Wigner's little group element W. The standard reference 3-momentum is directed along the 
z-axis of an arbitrarily chosen reference frame, with its x and y axes defining the two linear polarization vectors 
bf 2 , respectively. Hence kg = (1,0,0,1), and imposing the polarization gauge relates 3- and 4- polarization vectors, 

bi, 2 = (0,bi >2 ). 

The standard Lorentz transformation can be taken as 1 

L(k) = R(k)B z (k), (23) 

where B z (k) is a pure boost along the z-axis that takes kg to (k, 0, 0, k), and the standard rotation R(k) brings the 
z-axis to the desired direction k(#, <p) first by rotating by 9 around the y-axis and then by <j) around the z-axis, 

R(k) = R^RyiO). (24) 

The standard polarization vectors for an arbitrary momentum are defined as 

bi(k) = R(k)bf , (25) 

and a general real 4- vector of polarization can be written as 

f = cosxbi(k) +sinxb 2 (k). (26) 

Helicity is invariant under Lorentz transformations. The corresponding polarization vectors are 

b±(k) = -L(bx(k)Tib 2 (k)). (27) 

Under a Lorentz transformation A a state of a definite helicity acquires a phase e ±l ^ A,k ), which can be read-off from 
Wigner's little group element 

W(A, k) = L-\Ak)AL(k) = R z (0T(a,f3). (28) 

The little group element W leaves the standard momentum kg invariant. Here T{a, /3) form a subgroup which is 
isomorphic to the translations of a Euclidean plane and R Z {C) is a rotation around the origin of that plane, which in 
this case is also a rotation around the z-axis. 

If the transformation in question is a pure rotation 31, then the little group element is then a rotation i? z (£), and 
the polarization three- vector is rotated by 3? itself [30| • The phase x has a simple geometric interpretation: it is the 
angle between, say, bi(DUc) and 3ibi(k). It is zero if the standard polarizations of the new momentum are the same 
as the rotated standard polarizations of the old momentum. 

This is what happens if a rotation i? 2 (w) is performed around the current b 2 (k): the resulting phase Ax(i? 2 (w)) 
is zero [32J. Indeed, if k = k(0, 0), then by setting k' = i? 2 (w)k and using the decomposition 

R 2 {lo) = i?(k)i? y (w)i?- 1 (k), (29) 

we find that k' = k (6* + u>, <fi), and the little group element is 

W = R- 1 (k')R 2 (uj)R(k) = R~ 1 (8 + uj)R y (d)R y {oj) = 1, (30) 

indicating the absence of rotation with respect to standard polarization basis. 

In a curved space-time one has to provide the standard (xyz) directions at every point. The next section deals with 
this problem. 



1 When it does not lead to confusion we use the same letter to label a four-dimensional object and its three-dimensional part. In particular, 
R stands both for a Lorentz transformation which is a pure rotation, and for the corresponding three-dimensional rotation matrix itself. 
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B. Newton gauge 



The three-dimensional propagation equations (|19[) in stationary space-times result in a joint rotation of polarization 
and unit tangent vectors 

^ = f2xk, ^ = fixf. (31) 

In the Schwarzschild space-time test particles move in the plane passing through the origin 0, HIS^, O = —Eg x k 
and the covariant time-time component of the metric is related to the Newtoni an g ravitational potential ip as h = 
Soo = 1 + 2<p(r) — 1 — 2M/r. Requiring the resulting phase to be zero, as in (l8l. Il9l [26|. constrains a choice of local 
reference frames. We use the zero phase condition to make a physically motivated choice of standard polarizations 
that does not require references to a parallel transport or communication between the observers. 

Orient the local z-axis along the direction of the free-fall acceleration, z|| Vy> and the standard polarizations as 

y = b 2 := w x k/|w x k|, bi := k x b 2 , (32) 

where w is a free-fall acceleration in the frame of a static observer. Then f2 = Ob 2 , and the propagation induces no 
phase. 

In a general static space-time we take the z-axis along the local free fall direction as seen by a static observer. In 



the Kerr space-time its components are [34 



M(p 2 -2r 2 )\/A M M 2 ml 4n Mra 2 sin26> , sin 20 ml ^ , , 

^ = W-2Mr) ='^~- + W * = p3 (p2 , 2Mr) - + 0(r-% W , = 0. (33) 

This convention, that we will call the Newton gauge, is consistent: if we set z = — f in the flat space-time, then no 
phase is accrued as a result of the propagation. In addition to being defined by local operations, the Newton gauge has 
two further advantages. First, it does not rely on a weak field approximation to define the reference direction. Second, 
if the trajectory is closed or self-intersecting, the reference direction z is the same at the points of the intersection. 

C. Kerr space-time examples 

We impose the temporal gauge in local frames, hence 

a , 2Marf^ sin 2 (9 2Ma/ 3 sin0 , „ N 

f = Q^f = J — = L +0(r-% 34 

p z — 2Mr r z 

where we used / = f^e®^. This transversality reduces to a familiar three-dimensional expression k ■ f = k • f = 
f l k l + f 2 k 2 + fk 3 = 0. This also implies 

ft = ^)tF = e ( o)t/ 6 = 0. (35) 
In the Kerr space-time the Walker-Penrose quantity [4l| K 2 + iK\ , 

Ki =rB-aA cos 6, K 2 = rA + aB cos 0, (36) 

where 

A =(fc*/ r - k r f) + a{k r f' t > - k^f) sin 2 0, (37) 
B =((r 2 + a 2 )(fcV - k 6 ^) - a(fc 4 / 8 - k e /*)) sin0, (38) 

is conserved along null geodesies @ . Transversality and the gauge (|3"4")1 make the Walker-Penrose constants functions 
of only two polarization components, for example f 2 and f 3 . 

For a generic outgoing null geodesic in the asymptotic regime the constants become 

I<1 = Tout /out - ^out/out, K 2 = S/?out/out + Tout /out! ( 39 ) 
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where 

7out •— D esc O ut - a sin 6» ou t , (40) 
and the analogous expression gives the constant in terms of the initial data. Hence, 

/out = 0) /out = —52 T 2 (-g/^out^l + Tout ^2), /o 2 ut = ~ fl2 ! 2 ("^out^ - 7out-K"l)- (41) 

Pout i 7out Pout i 7out 

To determine the polarization rotation x we take the initial polarization to be, say, fi n = b 1 ", which often consid- 
erably simplifies the expression for K\ and A'2. Since K\ and Ki are linear functions of polarization, expressing the 
above result in the basis (b° ut ,t>2 Ut ) the desired rotation angle \ is observed from 

n ^f;Wr*Y (42) 



) \ ) V sinx 

where T is an orthogonal matrix that is described below. 

We assume fixed D and r\. Using Eq. (|32[) we find that at the limit n n — >-oo the initial standard linear polarization 
directions are 

bj n = 1 (0,sAnSinfl in ,£>), bf= — 1 (0,A-gAnSing in ). (43) 

'D 2 + A 2 sin 2 in ,/D 2 + A 2 sin 2 in 



Similarly, the final standard polarizations are 

b° ut = 1 = (0, s/3 out sin0 out , -£>) , fj° ut = 1 = (0, £>, s/3 out sin0 out ) . (44) 

^ 2 + /3 2 ut sin 2 in + /3 2 ut sin 2 # in 



The rotation matrices A^ n and iV ou t from the polarization bases to ((9, cf>) basis are given in the Appendix C. 

Following [2l1 ] using the asymptotic relationship between Walker-Penrosc constants and the components of polar- 
ization one can introduce a transformation matrix R, 

'40 (45) 



PI \ f? 

J out / \ ./ ii 



The transformation matrix has the form 



where the parameter x is given by 



^TTTfU-i 1 ' (16) 



AnTout - /?out7i] 



(47) 



with 7i n := Dcsc6 m ~ asin#j n . Finally, 
Expansion in the inverse powers of A results in 



7in7out + An/3 OU t 

T = N out RNr\ (48) 



T 21 =sin X = -^cosfl in + 0(A- 3 ). (49) 

This is our main new result. Few special cases are of interest. If the initial propagation is parallel to the z-axis and 
the impact parameter equals b, then 

A 2 = b 2 - a 2 , D = 0, (50) 

and the polarization is rotated by 

sin X=^2- + ^^+0(A- 4 ), (51) 
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and the antiparallcl initial direction gives the opposite sign. 

Motion in the equatorial plane corresponds to rj = (sin a = 0) and is qualitatively similar to the motion in 
Schwarzschild space-time. If the trajectory starts there but eventually moves outside, then the polarization is rotated 
by 

8M 2 a mli A , , , 

sinx = s— ^-sina + 0(A" 4 ). (52) 

In the case of initial propagation along the z-axis agrees with both [l8[ and [2(|, if we take into account the 
respective definitions of reference frames. On the other hand, in a generic setting our Eq. (|49[) differs by a power of A 
from 21 1 or (l9| (both works predict \ ~ A -3 but disagree on the pre- factors). The origin of this disagreement is in 



the following. In addition to using a more transparent series expansion, we use a different polarization basis from the 
one in [2l|. In particular, we do not require the knowledge of </> ou t — </>m- Moreover, [l|| considered only a Machian 
effect (Sec. HVj) . which is dominated by the reference- frame term [32j. 

IV. GEOMETRIC PHASE 

A more geometric perspective on polarization rotation is possible both in a static space-time, where we use the 
1+3 formalism, as well as in arbitrary space-times. While the phase A% is gauge-dependent for an open trajectory, 
we will show that it is gauge-invariant on a closed path. 

Consider the differential equations for polarization rotation J32|. By setting f = bi at the starting point and using 
the parallel transport equations Eqs. ([l} and ([2]) we arrive to the equation 

1 1 

-j? = V k (f.b 2 ) = T—f-V k b 2 . (53) 

dX cosx f-bi 

In a static space-time projection of the polarization on E3 results in 

d x _ 1 D(f-ba) 



d\ cos x dX 



(54) 



hence the desired equation is 




1 n Db 2 , . 

= wk+,^f — -. 55 
fb 2 dX 

The first term corresponds to the original Machian effect that was postulated in [23j], but the observable quantity 
involves both the Machian and the reference- frame terms (32| . 

Now we discuss the geometric meaning of these equations. First consider a basis of 1-forms (a 1 , a 2 , a 3 ) that is dual 
to the orthonormal polarization basis (bi,b 2 ,k) at every point of the tra ject ory. A matrix of connection 1-forms uj 
is written with the help of Ricci rotation coefficients uj 1 - as ui l - = ui\-a^ (El l39j . and a linear polarization is written as 

f = j^bi + / 2 b 2 = cosxbi + sinxb 2 . (56) 



Its covariant derivative equals to 



Df - 

— = ^?/ c b a+w ?-/ c k + k(/ c )b c , 0,0=1,2. (57) 



dX 3£J 3£ 
Taking into account the antisymmetry of the connections uj 1 - = uj^ = — u>^~, we find 

^ = (-br/ 2 + h 2 f l ) - u^kj + kfc(4/i + 4/ 2 ). (58) 

A comparison with Eq. (|3"T|) leads to the identification 51 1 = f2-bi = kui^, f2 2 = — and to an alternative equation 
for the polarization rotation, 



(59) 
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Given a trajectory with a tangent vector k one can define a SO (2) line bundle with the connection u> = ujl^kdX, 

similarly to the usual treatment of geometric phase {3^,|4(|. Freedom of choosing the polarization frame (bi,h>2) at 
every point of the trajectory is represented by a SO(2) rotation R^(ip(A)^J . Under its action the connection transforms 
asw^ RujR- 1 + R- 1 dR so 

_^u,.k + ^ + -. (60) 

In a static space-time we can consider a closed trajectory in space. Then the resulting phase is gauge-invariant, 
since the last term above is a total differential and drops out upon the integration on a closed contour, 

Ax = <j> uj-kdX + <j>u). (61) 

We can formalize it with the help of a curvature 2-form is introduced as 8 := dui + Hi A Hi. For the SO (2) bundle it 
reduces to 9 = dui, so by using Stokes' theorem the reference-frame term can be rewritten as a surface integral of the 
bundle curvature as 

(62) 

A more practical expression follows from our previous discussion: Ax = arcsinf OU f b2. Conservation of K\ and K% 
in the Kerr space-time ensures that if a trajectory is closed as a result of the initial conditions, then f out = fi n and 
Ax = 0. The Newton gauge is designed to give a zero phase along any trajectory in the Schwarzschild space-time. 
As a result of the gauge invariance of Eq. (|6"Tj) , no gravitationally-induced phase is accrued along a closed trajectory 
in the Schwarzschild space-time, regardless of the gauge convention. 

In a general space-time we introduce an orthonormal tetrad such that k = keo + fce3 at every point of the trajectory. 
We again impose a temporal gauge and set ei 2 = bi^, where the local polarization basis is chosen according to some 
procedure. Then from Eq. (f2"6"| it follows that 



V k f = k^4,f + f^e w = (-/ 2 e (1) + f\ 2) ) - (^ + 4)fc 

+ + 4)/ £ ) e (o) + + 4)/ £ ) e (3) - c =^ 2 - (63) 

where we used again the antisymmetry of the connection. From the parallel transport condition Vkf = 0, we can see 
that 

(4+4)/ £ = (4+4)/ £ = °- ( 64 ) 

and 

d X = (^52 + uj^kdX := (65) 

We cannot have a closed trajectory in chronologically-protected space-time, so to obtain a gauge-invariant result 
we consider two future-directed trajectories that begin and end in the same space-time points. This layout is similar 
to the two arms of a Mach-Zender interferometer [4j. We also align the initial and final propagation directions of the 
beams and use the same rules to define the standard polarizations. Similarly to the previous case the curvature is 
just O = dCl, and the Stokes theorem gives the phase as 

Ax = / - / 0= jn= Q. (66) 



V. CONCLUSIONS AND OUTLOOK 



Rotation one ascribes to polarization as light propagates on a curved background depends on the gauge conventions 
that are used along the way. However, the closed-loop phase is gauge-independent (assuming we use the same 
reference frame for transmission and detection). Similarly to other instances of a geometric phase, the phase Ax 
is given by the integral of the (bundle) curvature over the surface that is bounded by the trajectory. The Newton 
gauge provides a convenient local definition of the standard polarizations. It is motivated by its relative simplicity 
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and path-independence of the reference frames. It is also the gauge in which the statement of a zero accrued phase 
along an arbitrary path in the Schwarzschild space-time is correct. 

Our results characterize the behavior of polarization qubits on curved backgrounds. Specific regimes and realistic 
scenarios, both on the Kerr background an beyond, are to be investigated farther. Of special interest are those 



scenarios that will be experimentally feasible in the near future, such as sending photons between satellites [171 ]. Once 
the expression for Ax is obtained, it can be represent it as the action of a quantum gate [13], similarly to the the 
special-relativistic scenarios (l4ll43]. This will allow us to use the full toolbox of quantum optics to design experiments. 
Understanding of the polarization phase in GR opens new possibilities for optics-based precision measurements [33 ]. 
both classical and quantum. 
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Appendix A: Some aspects of the Kerr space-time 
1. Chronometric tetrad 

We fix an orthonormal tetrad at every space-time point (outside the static limit) by demanding that eo is the 
four- velocity of a static observer, and the vectors ex and e 2 are proportional to the tangent vectors d r and dg. Their 
covariant components e( Q ) M are 



/ / 2Mr 2Marsm 2 6 \ ,„ s 

""> = {-f~ M -W^m) (A1) 



and 



e (1) = (0, P /VK,0,0), e (2) = (0,0,^,0), e (3) = ^0, 0, 0, psm6^ 
Setting E = — ko = 1 leads to 



A - a 2 sin 2 (, 



(A2) 



p 2 ~2Mrak'f'sm 2 e 

k = n , A3 

p 2 - 2Mr 



and k r is expressed from the null condition k 2 = 0, 



_ A[l - (fc e ) 2 (A - a 2 sin 2 6) - (k^) 2 sin 2 OA] 



(fc'r = ^ — a — 2 • hr • (A4) 

A — a z sin a 

The spatial components of the momentum in the coordinate and tetrad bases are related to k m = e Jk 1 * as 

k 1 = k r p/VA, (A5) 

k 2 = k e Pl (A6) 



while k 2 = k 2 = k m k,r, satisfies 



fc a = **psinfl A /- ^-5-, (A7) 

A — sin a 



,9 2Afr 2M ml , K , 



p 2 - 2Mr 
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2. Constants of motion 

The re-scaled (£7 = 1) z-component of the angular momentum is 

D = o \» r ( fc0 (° 2 + r2 )(p 2 - 2Mr ) - 2Afar(l - ^asin 2 ^)). (A9) 
p z — 2Mr 

The re-scaled Carter's constant 

77 := K - {D -a) 2 , (A10) 

where the constant K is most conveniently expressed as 

K := {a sin 6 -D/ sin 9) 2 + (p 2 k e ) 2 . (All) 

The asymptotic expressions when r — > oo and the momentum components arc fixed arc 

D = k^r sin 9 + 0{r- 1 ), (A12) 
r] = (k? cos 2 9 + k?) r 2 + a cos 2 9 (a(kl + k? - 1) - AMk % smd) + 0( r - 1 ), (A13) 

A = r^q + q+0{r~ 1 ), (A14) 
CO sa = -^L+0(,- 2 ). (A15) 



Scattering with the impact parameter b that measures the coordinate distance r cos 9 from the z-axis is most 
conveniently described with the help of a fiducial Cartesian system. Its axes are "parallel" to the fictitious Cartesian 
axes of the Bycr-Lindquist coordinates, and in the asymptotic region it is just parallel to the global Cartesian grid. 
We introduce the momentum components pi that are related to the spherical components ki the the usual relations. 
In this case for the initial momentum parallel to the z-axis we have 



L> = 0(r 1 - 3 ), A= 0> 2 -a 2 + O(?T 2 )- ( Al 6) 

Some care is needed in treating R(r) as a function of the initial conditions when n — > oo. Introducing the constants 
Co, ci, C2, we write it as 



where 



R( r ) = : r 4 - ci.r(r A + cyr{r + c rf, (A17) 



cj := (A 2 - a 2 )/r\ = (kf) 2 + (kf) 2 + 0(r~ 2 ) (A18) 
ci := 2M(A 2 + a 2 - 2Aacosa)/ri = 2M{c\ - 2afc i^ &m 01 ) + 0( r f 2 ), (A19) 
c := -a 2 A 2 sin 2 a = -a 2 ((/c 2 ) 2 cos 2 ^ + {k\) 2 ) + 0{r^ 2 ) (A20) 



Appendix B: Trajectories 
1. The r integral 



The integral over r in Eq. ([5]) is split into three parts, 



where 



r r dr 

31 := / = = 3^oo — 3?i — 3^2, (Bl) 



.7=r- (B2) 



12 



We present a corrected expression for Jii and evaluate the term Jioo to the fourth order in A. Starting from 3?oo in 
using the method of [25j we factor out (r — r m ; n ) and substite r = r min /.T, writing R(r) as 

*■ = ^ fcin + ^ 2 + ^ + =: %^(*), (B3) 

■ L \ 1 min ' m in / 

where the constants A, B, C arc determined by the polynomial division. Noting that i?(l) = we have 

r 2 min = -A-— (B4) 

''min ^min 



which gives 



, , B C A 7 B Cx 
R(x) = -A + Ax 2 + + B5 

' mm ' mm ' mm ' m i n 



Rewriting the equation and inserting the values for A, B and C gives us 

, \ r 2 ■ , . 9w 9s / a 2 r? , 2M (D — a) 2 + n 1 + x + x 2 , 
R(r) = -^-(A - a 2 1 - x 2 1 + 2 . ' 5-1 + x 2 ) A , , T B6 

^ 4 V ^m( A -« ) r min A 2 - fl 2 1 + X 

Defining 



we finally get 



??(l+x 2 ) n (D^) 2 + f)l + i + I 2 

* W := 7a2 := -2 — 5 — , B7 

(A z — a z ) A^ — ar 1 + a; 



_ 1 

4= = ^=J== /'dx ((1 - x 2 ) (l + ^-f{x) + ^-9{x)\ \ 2 . (B8) 



lr miD VR VA 2 - a 2 7o 
Expansion the integrand up to the third order in A gives Eq. (15). 



1 / AM n /15, r2 3 2 . 2 \ 8aM 
•Anc = - 7T H h — rr — M a sm a -r- cos a 



M 
+ A3 



VA 2 - a 2 V A A 2 V 4 4 J A 2 

t-0(A" 5 ). (B9) 



128 

M 2 + a 2 (10 cos 2 a — 6sin 2 a) + 157rMacoso! 

3 



The second and third integrals can be performed as an approximation in powers of r\ , using Eq. (|A17[) . Since 

(c + ci_r)r\ \ 
r 2 (r 2 — c 2 r 2 ) / ' 



fl(r) = (r 4 - c 2 r 2 r 2 ) ( 1 + ^ t ^Ti ) , (BIO) 



the leading terms in the expansion are 

dr 1 C2ri ci / 2 2r. 2 — c|r? 

■ arcsm - 



r, c 2 ri r< 2 ^(c 2 ^) 2 c^r.r^r 2 -^ 

with ^ being either n or r 2 , where in the latter case we assume that ^2 > r%. 



+ o(n b ), (Bii) 



2. The integral 

We calculate the integral J l/\/~QdO using the method of @. Noting that 

6 = T] + (a-D) 2 - (asind - Dcsc6) 2 , (B12) 
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we perform a change of variables /i := cos 9 and obtain 



where fro r\ > 

e„ = aV_ + M 2 )04 - /i 2 ), < < M+: (B14) 

and 

/4 := ^ (V(A 2 - a 2 ) 2 + 4^ T (A 2 - a 2 )) . (B15) 

The calculation is performed using the approximations to the auxiliary integral 

r "+ dp 1 



9 p aJn\+n 2 _ 



F(iP,k 2 ), (B16) 



where F(ip, fc 2 ) is the elliptic integral of the first kind, and 

k 2 := /4/GU- cos^ := ^*/^+, (B17) 



with < tp < 7r and /i+ = + y 

The asymptotic expansion in the powers of A gives 



so 



/4 = sin 2 a + — sin 2 2a + C(A" 4 ) (B18) 
M_ = - cos 2 a+ -^prsin 2 2a + 0(A- 4 ) (B19) 



, o a 2 sin 2 a „ , . A , ,„ „, 

= ro — +OA" 4 , B20) 

A^ 



1 1 



1_ — S in 2 a + 0(A- 4 ) . (B21) 



Trajectories of the type 6\ 6* m j n /(9 max —> 62 correspond to cos# m i n = //+, and cos6* max = < 0, respectively. 
In both cases the integration leads to 

fl= , 1 --[F^ u k 2 )+F{^k 2 )\, (B22) 

and the expansion in powers of A leads to Eq. (fTU|) . 

Using that at the zeroth order (the flat space-time) + i\>2 = ""j we expand ^2 = t — ipi + as 

2 =tt-^i+^||. (B23) 

fc 

Equating the two expression for 31^ in the scattering scenario leads to 

6 = AM, (B24) 
£ 2 = 15Af 2 7r/4- 8Macosa, (B25) 
£ 3 = M(-6a 2 - 128M 2 +45Ma7rcosa- 24a 2 cos2a + a 2 (3 + cos 2^i) sin 2 a). (B26) 

Similarly, setting 62 = it — 0% + Sg, 

S * = T,% (B27) 

fc 
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and expanding the both sides of 



we find 



cos(# — Se) = /i+ cos ( arccos 1 — 5^ ] , (B28) 



*i = 4M Vsin " Q - C0S ^ (B29) 
smt^i 

In the special case of scattering with the initial momentum parallel (anti-parallel) to the z-axis the angular mo- 
mentum is zero, so (i+ = 1 and 9 nl = tpi B = 0, x, and it is easy to see that 9 ou t = 6^,, x — 6^, respectively. 

3. The (f> motion 

The equation for is given by Q 



R J aVr J A\Tr~ 

\J sm 2 eVe J VeJ 

If 9 const the first and the last terms cancel thanks to Eq. ((5J- For brevity we write the remaining terms as 

<f> 2 - fa =: ± J R <t '(r)dr± J T(9)d9, (B31) 
We again decompose the radial integral as 31* = 3^ — 3lf — 31* , where 

P OO /"OO 

31* = 2 / R*(r)dr, 3lf = / R*{r)dr. (B32) 



Following the same procedure as in Sec. B.l we find 



, _ a ( 8M-axcosa M(3Mx - 8a cos a) ,._ 3 . \ 
00 ~ JW~a^ I 4A + 6A2 +U[A > I ■ 



Expanding i?*(r) for r 3> fi we obtain 



*M = . .. ^ . . + (B34) 



a 2 D 



rvw^W) ( V 1 + + p^SD (r2 + ° 2 - 2Mr) 

C(r- 5 ), (B35) 













r 



r 2 \Jr 2 — cfrf 

where the coefficients Cj are defined in Appendix A2. Hence 



>(r)dr = - 2 ° M 2 2 j ' ( ^r V ~ c ^i + *f ~ r\ arctan ( g ) ) +G(r r 5 ). (B36) 



Integral J Td9 leads to the elliptic integral of the third kind. A standard change of variables cos 9 = \i = fi + cos t 
and the identity 

sin 2 6> = (1 - /4)(i +psin 2 t) (B37) 
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where 



1-/4 V A 2 

lead to 



p := -^2- = tan 2 a I 1 + ^- ) + C(A -4 ) (B38) 



T(6)d6=^^ / * d * (B39) 

aJuX+n^-P+Jo (l+psin 2 t)Vl-A: 2 sin 2 t 



with cos-0* = cos 6**/^+ as in Appendix B2. Hence 

A cos a 1 



T(0)d0 = _^ri^^_n(-p,^,fc 2 ), (B40) 



where ip, k 2 ) is the elliptic integral of the third kind. Expanding the pre-factor in the powers of A gives 

D 1 A sec a 



+ G(A- 4 ), (B41) 



while 

1 /7T 



and 



n(-p,V,0) = ^_ (|-arccot(yi+^tan^) =: Do(p,^) (B42) 



n(-p, fc 2 ) = n (p, v) + (V - n (p, V)) + <^ 4 )- (B43) 

2p 



Taking into account that 

1 



VT + P 1 V A 2 

we get for the scattering scenario 

D ^ $ 



cosa| l--2sin 2 a +0(A~ 4 ), (B44) 



where the first two terms of the series are 



out-0in = 7r— a4> (B45) 



A. 

cos 2 ?/' cos a + sin 2 ?/' sec a ' 



'!'. = ... 9 ~ _ ; , 2 , , 0*46) 



and 



£2(3 + cos2a) sec 3 a + 4seca(£? sm2V>i — £2 cos2?/>i) tan 2 a , 2 

$2 = - 7; o 7T< 1 - J a cosa + 2Ma, (B47) 

(1 + sec^a — cos 2tp tan z a)^ 



where £j are given in Appendix B2. 



Appendix C: Newton gauge relationships 



In this Appendix we present a relationship between the components of polarization in the chronometric tetrad 
and the Newton gauge basis. This relationship allows to calculate the Walkcr-Penrose conserved quantity from the 
operationally meaningful polarization information and in the asymptotic regime gives the matrices N of Eq. (|48[) . 
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The transversally of polarization k ■ f = ensures the linearity of the relationship between (/ 2 , / 3 ) and (f x , f y ), 
where 

f / »>, • /"!>,. (CI) 
Since h y = w x k/|w x k| and h x = — b y x k, by using Eq. ([55)1 we find 

by = jj. (w g fc §! -w 2 k % , -w* 2 ki + wjfcg) , (C2) 

^ = J\fk (^^1^2 — + , it;2(fc| + ^j) — wjfcj/c^, — ^3(^1^1 + ^S^g)) i (03) 

AA 2 = (u^) 2 + (u^) 2 + (-iwjfej + ajfc 2 ) 2 . (C4) 

Projecting f on these directions gives 



(C5) 

r = TTTT ('Wi - ^1*2) + / 3 (-^s)) ■ (C6) 



AO; 1 

We define the transformation matrix N as 

[f'J v/ 

In the asymptotic regime where y 1 — ^ it becomes orthogonal. Taking k m = (1, s/3/r, D/r sin#) we obtain in the 
leading order 1/r 

>+^ 2 /sin 2 V -0 D /^ 9 ) 
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